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Abstract - Voronoi cells of a discrete set in Euclidean space are known as gener- 
alized polyhedra. We identify polyhedral cells of a discrete set through a direction 
cone. For an arbitrary set we distinguish polyhedral from non-polyhedral cells us- 
ing inversion at a sphere and a theorem of semi-infinite linear programming. 
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1 Introduction 

The Voronoi diagram of a finite set in the n-dimensional Euclidean space 
is a popular concept in Discrete and Computational Geometry, cf. Aichholzer 
and Aurenhammer [ ] or Okabe et al. [ ], as well as in Minkowski Geometry, 
cf. Section 4 in Martini and Swanepoel [ ]. 

A natural generalization from a finite set is the concept of a discrete set. 
By definition, a subset of is discrete if its intersection with any bounded 
set of E^ is finite. Here M C E^ is bounded if sup^^^^^ ||x — y|| < oc 
with the Euclidean norm || • || based on the Euclidean scalar product ( • , • ). 
Equivalently, a subset of E^ is discrete if it has no accumulation point. 

We study the cardinality of half spaces needed to describe a Voronoi cell. 
A closed half space is defined for non-zero G E^ and A G M by 

H-{u,X) := {xgE^ I {x,u)<X}. 

The Voronoi diagram of a non-empty generator P C E^ is the tessellation 
of E^ consisting of the Voronoi cells 

V{p) := {x e¥P I \\x - p\\ < \\x - q\\ for all q eV}, p eV. 

By translational invariance we assume in this article that the origin 0^ of 
E^ belongs to V and we restrict to the cell V := V{On) at 0^^. Notice the 
closed half space representation 

2V = r]pev^{o„}H-{p,\\p\\'), (1) 
the intersection over the empty index set being understood as E^. 
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Figure 1: The Voronoi diagram of the generator V {{ — 1, z) \ z G Z}U{02} 
is depicted. The direction cone V = {(x, y) ^M? \ x < 0}U{02} is not closed, 
hence it is not finitely generated. Theorem 1.2 concludes that the Voronoi 
cell V is not a polyhedron. 



To catch the structure of a Voronoi cell of a discrete set we use the 
following definitions. Let a subset M C be given. The affine hull aff(M) 
resp. positive hull pos(M) of M consists of sums AiXi + • • • + Xn^n such 
that for i = 1, . . . , we have G M, G M and Ai + • • • + Xn = 1 resp. 
Xi > 0. Notice afr(0) = and pos(0) = {0^}. The convex hull of M is 
conv(M) := aff (M) H pos(M). If M = conv(M) then M is convex. If for 
any A > and AM := {Am | m G M} we have AM C M, then M is a cone. 
A convex cone is finitely generated if it is the positive hull of a finite set. A 
polyhedron is the intersection of finitely many closed half spaces. A bounded 
polyhedron is a polytope. If M is convex and if any intersection of M with 
a polytope is a polytope, then M is a generalized polyhedron. 

Remark 1.1. It is clear from (1) that a Voronoi cell of a finite generator is 
a polyhedron. It is proved by Gruber [3], Chapter 32, that a Voronoi cell of 
a discrete generator is a generalized polyhedron. 

The existence of a non-polyhedral Voronoi cell for a discrete generator is 
demonstrated in Figure 1. One of us has characterized a polyhedral Voronoi 
cell through the direction cone^ 

V := pos('P). 



Theorem 1.2 (Ina Voigt [ ]). IfV is discrete then the Voronoi cell V is a 
polyhedron if and only if the direction cone V is finitely generated. 

Further examples to apply this theorem are in Figure 2 and Figure 3. 
For subsequent discussions we recall that a convex cone is finitely generated 
if and only if it is a polyhedron, cf. Ziegler [ ]. If F is a convex subset of a 

"^In optimization the cone V is called the cone of feasible directions. 
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convex subset C C such that for x, y G C and < A < 1 the inclusion 
Ax + (1 — X)y G F always implies x^y ^ then F is a face of C. The zero 
dimensional faces are the extreme points and an extreme ray is a face which 
is a half line emanating from the origin. We notice that a polyhedron has 
at most finitely many faces, cf. §19 in Rockafellar [ ]. 

In this article we recover Theorem 1.2 from results about a generator V 
not necessarily discrete. As a tool we use the diffeomorphism of inversion 
at the unit sphere 

inv : \ {04 \ {04, x 

the reciprocal TZ iiw(V \ {^n}) and the convex reciprocal 

e := conv(7^U {0^}). 

Illustrations of C are given in the following figures. We apply in Section 2 
a theorem of semi-infinite linear programming by Goberna and Lopez and 
prove that the Voronoi cell V is a polyhedron if and only if the closure C is a 
polyhedron. The intuition is that V is completely surrounded by the Voronoi 
cells corresponding to the extreme points of C while V has no extension in 
unbounded directions of C. The surrounding cells are finite in number only 
if C is a polyhedron. In the case of a discrete generator V the condition 
relaxes to the condition that C is a polytope. 

Compared to the concept of convex reciprocal, the direction cone 

V = pos('P) = pos(7^) 

is nearer to the geometry of the generator V. We think of V as the area 
0^ + P occupied by V from point of view of O^?,. For arbitrary V we find in 
Section 3 that a polyhedral closure V is necessary for a polyhedral Voronoi 
cell V, i.e. a finite number of extreme rays of T> is necessary. On the other 
hand, a polyhedral V is not sufficient, a simple example being Figure 4. 

Stronger conditions apply to the case of a discrete generator we re- 
cover Theorem 1.2. The generator must be enclosed in a half space H~{u, 0) 
for non-zero u ^ to realize a non-poly heral Voronoi cell V. Then any 
asymptotic direction of the generator not in V, i.e. any accumulation point 
of I X G P \ {On}} not in P, makes a finite half space representation of 
V impossible, see Figure 1 and 2 for examples. (For the general generator 
there is no such condition, see Figure 5 as a counterexample.) 

In Section 4 we discuss the polar P* of the direction cone V. This is the 
normal cone at 0^ of the convex hull conv(P) of the generator, 

P* = { X G I (x, < for ah y G conv('P) }. 

As we noted above, the cone V = P** is useless as a sufficient condition for 
a polyhedral Voronoi cell. However, it is useful to decide if a Voronoi cell 
is bounded. This problem is resolved in the literature for finite generators, 
see e.g. Okabe et al. [ 
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2 The characteristic cone 



We apply a theorem from semi-infinite linear programming to the special 
case of the Voronoi cell V and obtain conditions on the cone of inequalities 
for the half spaces representation of V. The result is interpreted in terms 
of the convex reciprocal C. The cone of inequalities is well-known from the 
lifting construction for Delaunay triangulations. 

The starting point is the, possibly infinite, system of linear inequalities 

a {{p,x) < llpf \peV} 

satisfied by an unknown x E E^, if and only if x belongs to the doubly sized 
Voronoi cell 2V, see (1). The characteristic cone of a is^ 

K pos({(||;;|.) \per\{On}}u {{\-)}) C E"+i. 

The trivial equation for On ^ V is omitted. As a special case of Theorem 
5.13 in Goberna and Lopez [ ] the following equivalence holds. 

Theorem 2.1 (Goberna and Lopez). V is a polyhedron if and only if the 
closure K of the characteristic cone K is a polyhedron. 

In place of K C E^+^ we can study the convex reciprocal 

e = conv({x||x||-2 I x^V\{On}} U {0^}) C E^. 

We have K = |Ja>o ^i^^ {!}) C x {1} = {x^+i = 1}. We must be 
careful in a discussion of closures: if C is unbounded then K has non-zero 
points in the hyperplane Xn+i — but K does not. 

Proposition 2.2. The following statements are equivalent: 

(i) V is a polyhedron, 

(a) K is a polyhedron, 

(Hi) C is a polyhedron. 

Proof. We recall from Theorem 11.5 in Rockafellar [ ] that a closed convex 
set is the intersection of the closed half spaces that contain the set. More 
specific assertions about closures of convex hulls and closures of positive 
hulls are Corollary 11.5.1 and Corollary 11.7.2 in the same reference. As 
a consequence we can write for the same (possibly empty) index set / C 
E- \ {0,} X R 

^In [ ] the characteristic cone is used with the opposite sign (reflected at the origin) 
compared to our definition. 
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Figure 2: V := {|(cosh(x) + sinh(x), cosh(x) — sinh(x)) | x G ^Z} U {O2} 
is depicted with its Voronoi diagram (left). The Voronoi ceh V is not a 
polyhedron by Theorem 1.2 as P = {{x^y) G | x,^ > 0} U {O2} is 
not finitely generated. A second argument that V is not a polyhedron is 
Corollary 2.4 because the convex reciprocal C, having the infinitely many 
extreme points 7^, is not a polytope (right). 



and 

K = H-{{'_1),0) n r]iu,X)eiH-{{-\),0). 

Notice that the closed half spaces for the convex cone K have the origin 
On+i on their boundary hyperplane. A particular result is that K and C are 
both polyhedra or they are both not. Theorem 2.1 completes the proof. □ 

Proposition 2.2 is explained with non-discrete examples in Figure 4 and 
Figure 5. The statement of the proposition simplifies in the discrete case 
with the following remark. This is stated in Corollary 2.4 and explained 
with two discrete examples in Figure 2 and Figure 3. 

Remark 2.3. If P is discrete then the convex reciprocal C is compact and 
the characteristic cone K is closed. Observe that the reciprocal TZ is bounded 
having 0^ as the only possible accumulation point. Then TZ U {On} is com- 
pact and from Caratheodory's theorem follows that C = conv(7^ U {On}) is 
compact. Under^the hnear map a : E^^-*^ E^^-*^, (x, A) (Ax, A), the 
closed cylinder C := C x {A > 0} C E'^^^ is mapped to K. The kernel 
E^ X {0} of a does noi^ contain^ the direction y = (0^^, 1) of recession of C, 
i.e. a direction y with C + y C C, so K is closed, cf. Theorem 9.1 in [ ]. 

Corollary 2.4. IfV is discrete then the following statements are equivalent: 

(i) V is a polyhedron, 

(a) K is a polyhedron, 

(Hi) C is a polytope. 
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Figure 3: The Voronoi diagram of P := {(t, f{t)) \ t G |Z}U{02} is depicted 
(left). The graph of the real function f \ t ^ \{t + 2)^ has tangents through 
(—2,0) and (2,4) meeting O2. Thus, there is a steepest line through O2 that 
meets a point of {O2} so the direction cone V is finitely generated. Then 
Theorem 1.2 proves that V is a polyhedron. A second argument that V is a 
polyhedron is Corollary 2.4. The two branches inv(t, f{t)) are concave near 
O2 (for large |t|), so almost all points of TZ belong to the interior of C, which 
therefore is a polytope (right). 



Remark 2.5 (Delaunay triangulations). A Delaunay diagram of a finite 
generator V is defined as a tessellation of conv(7^) where a circumsphere 
of a cell is an empty sphere. Thereby a circumsphere of a cell is an empty 
sphere (also empty circle) if the interior of the corresponding ball has an 
empty intersection with V. A well known construction method for Delaunay 
diagrams is the lifting construction based on the map L : ^ E^^-*^, x ^ 
(x, ||x|p). A Delaunay diagram of V is obtained as the orthogonal projection 
of the lower convex hull of L{V) onto the (xi, . . . , x^?,) -plane. That is, the 
edges of 

{{x,z) G conv('P) X E I 3y < z such that {x,y) G conv(L('P))} C E''^^ 

projected orthogonally to the (xi, . . . , x^)-plane, produce a Delaunay dia- 
gram for 7^, see for example Okabe et al. ["']. In particular, the edges of 
pos(L(7^)) correspond to the edges emanating from 0^^ in this diagram. 

Heuristically, we consider the lifting construction for an infinite discrete 
generator V. The positive hull of the lifted generator is 

pos(L(P)) = pos((i--/-)), xg7^\{0,}) = U A[conv(7^)x{l}]. 

A>0 

Similarly as in Remark 2.3 we have pos(L(7^)) — ii V is unbounded 
(otherwise V is finite). If the Voronoi cell V is not a polyhedron then by 
Theorem 2.1 the characteristic cone K is not a polyhedron, so pos(L(P)) is 
not a polyhedron. This is in accordance with the Delaunay diagram having 
infinitely many edges emanating from the origin. 
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3 The direction cone 



We compare the convex reciprocal C to the direction cone V. While a 
polyhedral Voronoi cell V was found equivalent to a polyhedral closure C in 
the last section, we will see in this section that a polyhedral Voronoi cell V 
is a stronger condition compared to a polyhedral direction cone V. In the 
case of a discrete generator V these conditions are equivalent in accordance 
with Theorem 1.2. 

Let us study what consequences a polyhedral direction cone V can have 
for the Voronoi cell V. We continue in Remark 3.1 with a proof sketch of a 
necessary assertion, omitting to explain the concepts needed for a proof. 

Remark 3.1. If G C is a polyhedron containing the origin 0^ and 
if i7 = pos(G) then there exists an e > such that for the open ball 
S^(O^) := {x \ \\x\\ < e} we have 

G n BM = H n BM. (2) 

For a proof we can use the gauge of G defined for x G by 

7(x) := inf {A > I X G AC}. 

For (2) to hold it is sufficient to find some e > such that j{x) < holds 
for all X in the unit sphere S{H) := {x G | = 1}. The gauge 7 is a 
positively homogeneous function, whence it has finite values on H. On the 
other hand, with On G G, the positive hull H = pos(G) is a polyhedron and 
as such, is locally simplicial. These facts can be found in Rockafellar [ ]. 
Theorem 10.2 in this reference concludes that 7 is upper semi-continuous on 
H. Since the unit sphere S{H) is compact, 7 has a finite maximum there. 
This proves (2). 

Asking for sharpness of (2), let G := {x G E^ | ||x - (-^,0)|| < ^}, this 
is a closed disk touching the origin. Here the gauge 7 still is upper semi- 
continuous on. H \ {O2} = {{x,y) G E^|x < 0} but the unit sphere S{H) is 
not compact and for each e > we have GnB^{On) $ H nB^{On)- A related 
example is the convex reciprocal for the example in Figure 1 with reciprocal 
7^ = {(-1, t)(l + 1^)-! I t G Z} included in the boundary of G. 

Proposition 3.2. If the generator V is discrete and if the direction cone V 
is a polyhedron, then the convex reciprocal C is a polytope. 

Proof. If the cone V = pos('P) is finitely generated, then we can assume 
that it is finitely generated by points of the generator V or by points of the 
reciprocal 7^, likewise: for ri, . . . , r/^ G 7^ we have 

V = pos(ri, . . . ,r/e). 
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Figure 4: A finite sketch of the Voronoi diagram for unit sphere with center 

V := {(cos((/p), sin((/^)) | (/p G [0, 27r)} U {O2} is depicted (left). As the convex 
reciprocal C — C — conv(7^) is not a polyhedron, the Voronoi cell V is not 
a polyhedron, cf. Proposition 2.2. Another example is the Voronoi cell V of 

V := {(e~^, 1) I X = 1, 2, 3, . . .} U {O2, (0, 1)} (right). Here the convex recip- 
rocal C, having the infinitely many extreme points 7^, is not a polyhedron 
(middle). Still, the direction cones are polyhedral in both examples. 

Now we consider the polytope 

C := conv(0^,ri, . . . ,r/e). 

Since V — pos(C) we meet the assumptions of Remark 3.1 and can infer that 
D\inv(C) is bounded. So this set contains at most finitely many points of the 
discrete generator V. Then all but finitely many points 5i, . . . , 5/ of 7^ belong 
to C and therefore the convex reciprocal C — conv(072, ri, . . . , r/^^, 5i, . . . , 5/) 
is a polytope. □ 

Corollary 3.3. // the generator V is discrete and if the direction cone V is 
a polyhedron then the Voronoi cell V is a polyhedron. 

The above conclusion follows from Proposition 3.2 and Corollary 2.4, 
Figure 3 shows an application. Figure 4 demonstrates that a polyhedral 
direction cone is not sufficient for a polyhedral Voronoi cell without the 
assumption of the discrete generator. 

Assuming the discrete generator, a polyhedral closure V of the direction 
cone is not sufficient, we remember Figure 1 and Figure 2. Now we will see 
that a polyhedral closure V of the direction cone is necessary regardless of 
the generator. 

Proposition 3.4. // the convex reciprocal C is a polyhedron then the direc- 
tion cone V is a polyhedron. If the closure C is a polyhedron then the closure 

V is a polyhedron. 

Proof. We have the trivial chain of inclusions 

V = pos(C) C pos(C) C pos(C) = V. 



8 



Figure 5: We consider the generator V := {^2^Pa{i)^Pa{2)^Pa{?>)^ - - -} with 
Pa := 2cos(a)(— sin(a),cos(a)) and a{t) := |e~^ The closed convex recip- 
rocal C is a triangle (middle), so the Voronoi cell V is a polyhedron by Propo- 
sition 2.2 (left). While the direction cone V = pos({(A, 1) | -1 < A < 0}) is 
not closed, the closure P is a polyhedron by Corollary 3.5. The idea to the 
example (right): we narrow V below certain lines through (0, 1). 



With C being a polyhedron containing the origin. Corollary 19.7.1 in [ ] 
proves that the positive hull pos(C) is a polyhedron. Thus, with pos(C) be- 
ing closed we obtain that V — pos(C) is a polyhedron. Assuming that C is 
a polyhedron, we can replace C by C and obtain V = pos(C). □ 

Corollary 3.5. // the Voronoi cell V is a polyhedron then the closure V of 
the direction cone is a polyhedron. 

Corollary 3.6. // the generator V is discrete and if the Voronoi cell V is a 
polyhedron then the direction cone V is a polyhedron. 

The above conclusions follow from Proposition 3.4 together with Propo- 
sition 2.2 and Corollary 2.4 in this order. Figure 5 introduces a polyhedral 
Voronoi cell V where V is not closed. So, a polyhedral direction cone V 
is necessary for a polyhedral Voronoi cell V only in the case of a discrete 
generator. 

4 Bounded cells 

Using normal cones we prove a condition when the Voronoi cell V is bounded. 
For a discrete generator V this is a condition when V is a polytope. 
The normal cone of a convex subset C C at x G C is defined by 

N{C, x) := { G I {y-x,u) <0 for all y e C }. 

This is the set of vectors G E^ having no acute angle at x with any point 
y ^ C. An example is 

7V(conv(7^),0^) = N{C,On) = {u eE"" \ {y,u) < for all y e V}, (3) 
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where the set equahties hold because for G the inequahty {y^u) < 
for sll y ^ V is equivalent to this inequality for all y G 7^ or for all y in 
the convex hull of one of these sets. The inequalities are even equivalent to 
these with y running through the positive hull pos(V) = P, so the polar of 
the direction cone P* := { G E^ | i/) < for all ^ E P } satisfies 

V* = iV(conv(7'),0„) = N{C,On). (4) 

The recession cone of V describes unbounded directions of V, it is 

rec(V) := {ueE'' \ x + A^x G V for all x E V, A > }. 

Lemma 4.1. The equality of cones A^(conv(P), 0^^) = rec(V) holds. 

Proof. We have the representation 2V = PipG'PxiOn} ll^lP) ^^^^ 

spaces (1). With Corollary 8.3.3 in [6] the recession cone rec(V) = rec(2V) 
becomes 

r]pev\{On}^^<H-iP,\\p\\')) = r]pev\{On}H-iP,0). 
The proof is completed by (3). □ 

As a consequence we can determine boundedness of the cell V. For 
M C E^ let M° denote the interior of M in the topology of E^. 

Proposition 4.2. The Voronoi cell V is bounded if and only if On E C° if 
and only if On E conv(P)°. IfV is discrete, then V is a polytope if and only 
if one of these equivalent conditions holds. 

Proof. By Theorem 8.4 in [ ] the Voronoi cell V is bounded if and only if the 
recession cone rec(V) is zero. Using the equality A^(conv(7^), On) = rec(V) in 
Lemma 4.1 this is equivalent to a zero normal cone A^(conv('P), 0^^). Now by 
Theorem 13.1 in [6], the normal cone at 0^^ is zero if and only if 0^ E conv(P)° 
holds. Using (3) we can argue with C in place of conv(P). In the discrete 
case, the cell V is a generalized polyhedron, see Remark 1.1. But a general- 
ized polyhedron is bounded if and only if it is a polytope. □ 

The zero normal cone responsible for a bounded Voronoi cell in Propo- 
sition 4.2 is equivalent to the equality P = E^ through polarity of closed 
cones, cf. (4) and Theorem 14.1 in [ ]. Since a convex cone P C E^ is 
included in a half space, we have for arbitrary generator V 

V is bounded ^ P = E^ ^ P = E^. 

Figure 6 shows an example of a discrete generator in E^, which is bounded 
in ^/-direction and where nevertheless every Voronoi cell is a polytope. 
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Figure 6: All Voronoi cells of the generator V '.— 1 — ^) | n E N} are 

polytopes by Proposition 4.2 as conv(P) = {{x^y) G | \y\ < 1} is open. 



5 Conclusion 

By considering a Voronoi cell V as a problem in linear semi-infinite program- 
ming, we obtain an equivalent condition when V is a polyhedron. This is 
the condition that the closure C of the convex reciprocal is a polyhedron and 
the condition simplifies in the case of a discrete generator V to the condition 
that the convex reciprocal C is a polytope. 

While the polyhedral Voronoi cell V implies the polyhedral closure V of 
the direction cone, the polyhedral cone TD follows only in the discrete case. 
Conversely, the polyhedral direction cone V implies a polyhedral Voronoi 
cell only in the discrete case. 

The closure V gives only advice if the Voronoi cell V is bounded or not. 

Acknowledgment. Special thanks to Andreas Knauf for general support 
during our meeting in Erlangen in August '09, to Oliver Stein who brought 
[ ] to our attention and to the anonymous referee, who brought [ ] to our 
attention. 
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